Correlation effects on resonant tunneling in one-dimensional quantum wires 
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We study resonant tunneling in a Luttinger liquid with a double barrier enclosing a dot region. 
Within a microscopic model calculation the conductance G as a function of temperature T is deter- 
mined over several decades. We identify parameter regimes in which the peak value Gp{T) shows 
distinctive power-law behavior. For intermediate dot parameters Gp behaves in a non-universal way. 



A decade ago resonant tunneling through a double bar- 
rier embedded in a Luttinger liquid was intensively stud- 
ied theoretically. 1, ^] At temperature T = 0, in the 
limit of an infinite system, and for symmetric barriers 
the conductance G was predicted to show a series of in- 
finitely sharp resonances if the gate voltage Vg applied to 
the dot region is varied. [ij For spinless fermions, on which 
we focus, the peak conductance Gp at resonance voltage 
Vg is given by /h if the interacting wire (denoted "the 
wire" in the following) is connected by "perfect" contacts 
to non-interacting semi-infinite leads as will be the case 
here. Increasing T resonances are still present but acquire 
a non- vanishing width w{T) oc T^~^ , 1] with the inter- 
action dependent Luttinger liquid parameter K {K = 1 
for non-interacting particles). We here restrict ourselves 
to the case 1/2<K< 1. For AVg = \Vg-V^\ ^ 0, G(T) 
scales to zero as T^"^,||l| with ub = 1/K — 1. This also 
holds for an asymmetric double barrier regardless if the 
conductance is evaluated at or off resonance. 

Using second order perturbation theory in the barrier 
height the T > deviations of Gp{T) from /h were 
determined to scale as r^^.,2j At further increasing T 
a regime of uncorrelated sequential tunneling (UST) was 
predictedjj] based on a perturbative analysis in the in- 
verse barrier height, where Gp{T) oc T"^~^ and w{T) oc 
T. This regime is bounded from above by the level spac- 
ing of the dot Ad = ttvf/Nd, where vp denotes the 
Fermi velocity, Njj the number of lattice sites in the dot, 
and the lattice constant was chosen to be 1. Within the 
same perturbative approach and for Ad < T < B, where 
B is the bandwidth, Gp{T) increases as T^"^ for increas- 
ing In contradiction to the temperature depen- 
dence following from UST, in transport experiments on 
carbon nanotubes a suppression of Gp with decreasing 
T was observed for T ^ Aij.jBI Using another approxi- 
mation scheme "correlated sequential tunneling" (CST) 
was predicted to replace UST.[Q| It leads to Gp oc T-^"^^^ 
which was argued to be consistent with the experimen- 
tal data. This consideration has stimulated a number of 
theoretical works. 011 H ^ 

Quantum Monte Carlo (QMC) results'!^ for Gp(T) 
were interpreted to be consistent with CST for T ^ A^i. 

In Refs. a "leading-log" type of resummation for 

the effective transmission at the chemical potential was 
applied and led to important insights. In this approach 



no signature of CST was found, but it was criticized for 
only being valid in the limit of small interaction. That 
the "leading-log" method indeed does not capture all in- 
teresting interaction effects away from K \ can be seen 
considering the single barrier case for which the conduc- 
tance for different temperatures and barrier heights can 
be collapsed onto a single curve by using a one-parameter 
scaling ansatz.Tl While the resulting scaling function is 
known to depend on A', 1] the "leading-log" approach 
leads to the non-interacting K — 1 scaling function inde- 
pendently of the interaction strength chosen. 

In Ref. a A'^^i = 1 double barrier model in which 
chiral fermions carry the current was treated for K ~ 
1/2 by refermionization. After fine tuning the Coulomb 
coupling across the barriers Gp(T) was predicted to show 
the non-interacting behavior, i.e. e'^ /h — Gp{T) oc for 
T ^ and Gp{T) oc T^^ in the UST regime. This 
finding is surprising and not supported by any of the 
results obtained for more general models using the above 
mentioned methods. The physics of this particular model 
thus seems to be non-generic. 

We study the problem applying the fermionic func- 
tional renormalization group (fRG) method to the model 
of spinless fermions on a lattice with hopping t = 1 and 
nearest neighbor interaction [/ > at half filling, [l^ 

By using a method which can be applied for any bar- 
rier height and shape but is restricted to weak to interme- 
diate interactions our study is complementary to the one 
using perturbation theory in either the barrier heights 
or the inverse barrier heights. 0,3 Considering the sin- 
gle barrier case and comparing with exact asymptotic 
and numerical results we have shown previously that our 
method leads to meaningful results even for fairly large 
U £ 2.[il [H m In Ref. it was shown that for 
K — 1/2 our data for G collapse onto the exactly known 
scaling function of the local sine-Gordon model. With- 
out any further approximations the fRG can be applied 
to the double barrier problem and we thus expect that 
regimes with universal power-law scaling and the related 
exponents can reliably be determined. Accurate results 
can also be obtained for regions of dot parameters not in- 
vestigated so far in which Gp{T) shows non- universal be- 
havior, i.e. a behavior which cannot be characterized as a 
power-law scaling with an exponent expressible in terms 
of K. Studying these turns out to be relevant for the 
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interpretation of the QMC data|lO| and the experimen- 
tal results. (51 In contrast to the other methods which can 
only be applied in certain temperature ranges the fRG 
provides reliable results for G{T) on all energy scales. For 
large barriers and large as well as small dots we identify a 
variety of temperature regions in which Gp{T) shows dis- 
tinctive power-law behavior with different exponents ex- 
pressible in terms of K. For intermediate barrier height, 
intermediate dot sizes, and T ^ Ajj, Gp{T) decreases 
with decreasing T. The behavior of Gp{T) in this tem- 
perature regime depends on Nu and the height of the 
barriers and is thus non-universal. Whereas our results 
are consistent with the QMC data, we do not observe 
signatures of CST. 

The model is given by the Hamiltonian 

oc 
j = -oo 

with 1 -C < jr ^ N, in standard second-quantized 
notation. The lattice can be divided in three parts: (1) 
the wire with nearest neighbor interaction Uj across the 
bonds + with j € [1, (2) the dot (embedded 

in the wire) on sites j £ [ji,jr] {Njj = jr ~ ji + 1) which 
can be shifted in energy by the onsite energy Vg and is 
separated by barriers of strength Vj/^ from the rest of the 
lattice; (3) the non-interacting leads on sites j < 1 and 
j > N. Besides site impurities we also consider hopping 
impurities as barriers. In this case the interaction [in the 
last term of Eq. ^] and hopping [in the first term of Eq. 

across the bonds — ji and j,. , > + 1 are set to zero. 
Furthermore the second and third terms in the first line 
of Eq. are replaced by — t;c]^„;^c^^ — trCj^Cj^^-^ + h.c. 
with the reduced hopping matrix elements ti/r < 1. As in 
the other theoretical studies mentioned we avoid effects 
of the contacts. We thus model "perfect" contacts by 
a smooth spatial variation of the interaction Uj close to 
the sites 1 and N. Starting at site 1 in our calculations 
the interaction was turned on in an arctan-shape over the 
first 100 sites and similarly turned off approaching site 
A/'.[l^ For T = 0, vanishing barriers, and up to = 10^ 
this gives 1-G/ (e^//i) < 10^^. The constant bulk value 
of the interaction is denoted by U. 

In linear response and the absence of vertex corrections 
the conductance is given bv[lS| 

GiT) = -{e^/h) 1' |t(e,T)|2 dJ{e/T) de , (2) 

with \t{e,T)\^ = (4 - e2)|G'i,A,(e,T)|2 and / being the 
Fermi function. The one-particle Green function Gij has 
to be computed in the presence of interaction and leads. 
To this end we extend a recently developed fRG scheme, 



which has been shown to provide excellent results at 
T = 0,[ll[ll[l3 to finite temperatures. The starting 
point of this approach is an exact hierarchy of differen- 
tial flow equations for the self-energy and higher or- 
der vertex functions, with an energy cutoff A as the flow 
parameter. The hierarchy is truncated by neglecting n- 
particle vertices with n > 2, and the 2-particle vertex 
is parametrized approximately by a renormalized near- 
est neighbor interaction U^. The spatial dependence of 
the renormalized impurity potential, described by T,^, is 
however fully kept. This procedure is perturbative in the 
2-particle interaction, but non-perturbative in the impu- 
rity strength. Within our approximate scheme T,^ has no 
imaginary part, which implies the vanishing of the ver- 
tex corrections to the conductance formula (O-E^ The 
Green function Gi^n entering Eq. ^ is obtained from 
the Dyson equation with S^^°. 

As a test of our method at finite T we first discuss re- 
sults for the single barrier problem. We find the expected 
scaling behavior of the conductance (for an example see 
the solid line in the upper panel of Fig. I^J.^l, 12] From 
the suppression of G{T) oc T^"^ , the exponent as can 
be read off. Within our truncation scheme we obtain 
an approximation a§'^(t/) for asiU). Consistent with 
the exponents of the 1 /A^-scaling determined in Ref . 
at T = we find af^{0.b) = 0.165, ag^(l) = 0.35, 
and 0^*^(1.5) = 0.57 in good agreement with the ex- 
act exponents q;_b(0.5) = 0.1608, 0^(1) = 1/3, and 
Q!b(1.5) = 0.5398 known from the Bethe ansatz solution 
of the bulk model, [l^ 

For the double barrier problem we start the discussion 
with the case of strong, symmetric barriers. In Fig. ^ 
G{Vg) is shown for No = 6, U ^ 0.5, Vi/r = 10 and 
different T. In all figures the wire length is chosen to 
be iV = lO'^ sites in rough agreement with the length of 
carbon nanotubes accessible to transport experiments. Q 
For sufficiently small T resonance peaks are clearly devel- 
oped. For T — ^ 0, Gp{Vg) tends towards e'^/h. Because 
of iV < oo even at T = the peaks have finite width. For 
all peaks the T — width scales to zero as A^^"'^~"'^,0 
with the fRG approximation K^'^ for K. The different 
height and width of the peaks at fixed T is a finite band 
effect. In the following we always study the behavior of 
the resonance closest to Vg = 0. We have verified that 
Gp{T) for the other peaks shows a similar T-dependence. 

In the upper panel of Fig. El Gp{T) is shown for 
U = 0.5, Vi/r = 10, and No = 2, 6, 100. The lower panel 
contains the logarithmic derivative of Gp{T) from which 
the exponents of possible power-law behavior can be read 
off directly. Several temperature regimes can be distin- 
guished. For T larger than the bandwidth, Gp{T) oc T^^ . 
For T < Ad, marked by arrows in Fig.|3 down to a lower 
bound T* discussed below we find Gp{T) oc T"^ -i as 
indicated in the lower panel^f Fig.|21by the dashed line. 
This is the regime of UST, 4] in which the integral Eq. 
^ leading to Gp{T) is dominated by a single peak of 



3 



0.03 




FIG. 1: The conductance as a function of gate voltage for 
No=6,U = 0.5, V,/r = 10, N = W^, and different T. 
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FIG. 2: Upper panel: Gp(T) for (7 = 0.5, N = 10"*, Vi/r = 10, 
and JVd = 2 (circles), 6 (squares), and 100 (diamonds). The 
arrows indicate Ad. The solid curve shows G(T)/2 for a 
single barrier with V = 10 and U = 0.5, A'^ = 10*. Lower 
panel: Logarithmic derivative of Gp(T). Solid line: 2a§'^; 
dashed line: a^'^ — 1; dashed-dotted line: 2a§*^ — 1. 



|i(e,T)p, but T is still much larger than the width of 
this peak. Peaks in Gp{Vg) are clearly separated and 
the width w{T) oc T can be read off. The size of the 
UST regime grows with increasing No and shifts towards 
smaller T. For Nd > 4 we find non-monotonic behav- 
ior of Gp{T). For Nd ~ 30 the region of decreasing 
conductance (with decreasing T) can be described by a 
power-law with exponent 2a§^ (solid line in the lower 
panel of Fig.jJ).^^ For these T many peaks of \tie,T)\'^ 
contribute to the integral. For Vj/^ ^ 1 and Nd ^ 1, 
i.e. in the parameter regime with distinctive power-law 
scaling, the non-monotonic behavior can be understood 
using Kirchhoff's law. The solid line in the upper panel 
of Fig. 121 displays G{T)/2 obtained for a single barrier 
with V = IQ and U = 0.5. The comparison shows that 
the conductance of the double barrier problem for T up 
to the local minimum can be determined by adding up re- 
sistances of the related single barrier case. We do not ob- 
serve indications of CST with exponent 2a^ — 1, shown 
as the dashed-dotted line in the lower panel of Fig. |21 For 



T 0, Gp{T)/{e'^/h) saturates at 1. Due the finiteness 
of the wire for T <^ Aw, with Aw = nvp/N, we find 
1 - Gp{T)/{e^/h) oc T^. For small Nd and intermediate 
to small barriers at increasing T this behavior is followed 
by a regime with 1 - Gp{T)/{e^/h) oc T'^^'^° .2\ All the 
above scaling regimes we also find for U — 1 and 1.5. 

Temperature regimes with power-law behavior and the 
same scaling exponents were obtained using a field theo- 
retical model and perturbation theory in Vj/^ a-s well as 
l/^/r-@jl3l The agreement of the results obtained from 
two complementary approximate methods and for differ- 
ent models constitutes a strong argument for the correct- 
ness of the physics obtained. These considerations put 
very tight limits on the parameters for which a scaling 
regime with exponent 2q;b — 1 could occur. 6] As an ad- 
vantage of our method results for Gp on all energy scales 
down to T = can be obtained within one approxima- 
tion scheme and for more realistic lattice models. Fur- 
thermore the fRG allows for extensions to non-perfect 
contacts which have to be considered when comparing 
with experiments. 
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FIG. 3: Upper panel: Gp{T) for U = 1.5, iV = 10*, No = 10 
with Vi/r = 0.8 (circles), Nd = 30 with Vi/r ~ 1.5 (squares), 
and Nd = 50 with Vi/^ = 0.8 (diamonds). The arrows in- 
dicate Ad. Lower panel: Logarithmic derivative of Gp{T). 
Solid line: 2q|^; dashed line: q:|'^ - 1. 

We next discuss the case of weak to intermediate bar- 
riers, dots of a few ten lattice sites, and larger U. In this 
part of the parameter space we find non-universal behav- 
ior of Gp{T) which so far has not been studied theoret- 
ically and provides an alternative interpretation of the 
QMC data. ^3 It might also be relevant in connection 
with experiments since it is not clear that the experi- 
mental dot parameters fall in the regime of universal scal- 
ing. In Fig.iniGp(r) is shown for U = 1.5, Nd = 10 with 
Vijr = 0.8 (circles), Nd ~ 30 with V/y^ = 1-5 (squares), 
and Nd — 50 with VJ/^ — 0.8 (diamonds). These pa- 
rameter sets are chosen to be close to the set used in 
the QMC calculations, fioj In the limit of weak barriers 
and small dots (circles) neither the power-law with the 
exponent 2oiB nor the one with as — 1 are clearly de- 
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veloped. For increasing Njj or increasing Vi/r the UST 
regime gets more pronounced. For weak to intermedi- 
ate barriers no plateau with exponent 2a b is estabhshed 
(see the lower panel). The behavior in the regime of de- 
creasing Gp{T) (for decreasing T) depends on VJ/^ and 
Nu and is thus non-universal. It might be identified in- 
correctly as a power-law with an exponent significantly 
smaller than 2aB, i-e. in the vicinity of 2aB — 1- 

Our fRG results for Gp{T) in the case of asymmetric 
barriers and the off resonance scaling of G{T) for both 
symmetric and asymmetric barriers will be presented in 
Ref. ^3 • For later reference we mention that in all these 
cases we obtain T^"b'^ for T — > as expected. 0] 

Using the projection operator onto the subspace of the 
one-particle states of the dot and the orthogonal com- 
plement to the left and right of it the result |i(£,T)p = 
(4 — e^)|Gi^Ar(e, T)p can be written in a form which is 
useful for obtaining a deeper understanding of the dif- 
ferent scaling regimes discussed above and to work out 
the differences of our approach to the "leading-log" ap- 
proximation. As an example we take the case of hopping 
barriers and = 1. Then one obtains 

|t(£,T)p= (3) 

4r;(£,T)r,(£,r) 

[e-Vg- Qi (e, T) - ^rie, T)f + [Ti {e,T) + (e, T)f 

with T I = ImG° fli = ReG° and sim- 

ilar expressions for / — > r. The Green function G'^j{e,T) 
is obtained by considering as an effective potential 

outside the dot and setting ti/^ = 0. 

For Vg = 0, e -Vg - ni{e,T) - Vlr{e,T) vanishes at 
£ = and we find a resonance. Due to the generalized 
Breit-Wigner form of Eq. ^ for ti ~ tr the transmission 
at resonance is 1 independent of the dependence of F; 
on its arguments and the parameters. It is important to 
note that in this case and for T ^ 0, (i) F;/,r(£ — 0,T) 
does not follow the single barrier scaling T"^ but instead 
saturates (even for N oo). For asymmetric barriers 
with tr < U < 1 and at resonance Tr{£,T) oc 
while (ii) F;(£,r) increases as T~"b . Similar scaling 
is found for the e and 1/N dependence. These power- 
laws are cut-off by the largest of the scales T, e, and 
Aw- Inserting the expressions for Fj/^ into Eq. Q leads 
to a power-law suppression of |t(£,T)p with exponent 
2aB^ which directly gives Gp{T) oc T^"B mentioned 
above. We next study the off resonance case with Vg ^ 
and ti — tr < 1 including the single barrier situation 
{ti/r — 1). Because of the symmetry in the following we 
suppress the indices l/r. We focus on T = and consider 
Aw oc 1/N as the relevant energy scale. Then the inte- 
gral in Eq. Q can be performed and \t{0, 0)^ is directly 
related to G. For small Vg we (iii) find F(0,0) oc const, 
and Vg + 2f](0,0) oc VgiV^"^'^'^, which combines to the 
expected 1 - G/{e^/h) oc iv2(i-A:«'^)_ por large Vg, (iv) 



F(0, 0) oc t^N-"B and Vg + 2n{Q, 0) oc Vg which leads 
to G oc iV-2"s . The observations (i)-(iv) are not cap- 
tured by the parameterization of the effective transmis- 
sion |i(£,r)|2 used in the "leading-log" method. HHIIJ 
In this approach the important energy dependence of the 
real part f2 [see (iii) and (iv)] does not appear. 

Using Eq. Q and its generalization for Nu > 1 the 
UST regime can be understood quite simply. We here 
focus on symmetric barriers. Approaching from large T 
UST sets in at a scale at which only one peak of \t{e, T)p 
contributes to the integral in Eq. ||2J) (respectively for T 
smaller than the bandwidth at Nd — 1). As can be seen 
in Fig. 121 this scale is significantly smaller than A^i. For 
these T > |£| we find F(£,T) oc t'^T°'B'^ /No- Perform- 
ing the integral in Eq. |(2J) leads to Gp{T) oc T"b -i and 
w{T) oc T. This scaling continues down to the temper- 
ature T* oc (^V^d)^/^^""^'^^ at which the width of the 
peak in |i(£,T)p is equal to T. This defines the lower 
bound of the UST regime. 

In conclusion, applying the fRG we obtained a com- 
prehensive picture of resonant tunneling in a microscopic 
model for a Luttinger liquid with 1/2 < K <1. Depend- 
ing on the dot parameters we found distinctive power-law 
scaling of Gp (T) , but also regimes with non- universal be- 
havior. No signatures of CST were found. 
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